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ABSTRACT 

Archival Spitzer observations of 41 starburst galaxies that span a wide range in metallicity reveal for 
the first time a correlation between the [Fell] /[Nell] 26.0/12.8 /xm ratio and the electron gas density 
as traced by the 18.7/33.4 /jm [SIII] ratio, with the [Fell] /[Nell] ratio decreasing with increasing gas 
density. The correlations of the [Fell] /[Nell] ratio, the PAH peak to continuum strength & metallicity 
found in an earlier paper were confirmed for a larger sample. We also find a strong correlation between 
the gas density and the PAH peak to continuum strength. Using shock and photoionization models, 
we see the driver of the observed [Fell] / [Nell] ratios is metallicity. The majority of [Fell] emission in 
low metallicity galaxies may be shock-derived, whilst at high metallicity, the [Fell] emission may be 
instead dominated by contributions from Hn and in particular from dense PDR regions. However, 
the observed [Fell] /[Nell] ratios may instead be following a metallicity-abundance relationship, with 
iron being less depleted onto grains in low metallicity galaxies - a result that would have profound 
implications for the use of iron emission lines as unambiguous tracers of shocks. 
Subject headings: galaxies: starburst - galaxies: stellar content - ISM: lines and bands - infrared: 
galaxies 



1. INTRODUCTION 

The presence of massive stars within starbursts un- 
doubtedly plays a huge role in determining the physical 
conditions within the local interstellar medium (ISM). 
High-mass (Mi n it > 8M©) stars formed within typical 
starbursts drastically affect the dynamics of the sur- 
rounding ISM, through not only the release of ioniz- 
ing photons which destroy molecular material, but also 
via supernovae (SNe) which provide thermal and kinetic 
energy input into the ISM. The effects of photoioniza- 
tion by high-mass stars on the obs erved dearth of PAHs 
have recently been inve stigated dMadden et ail 120061: 
IWu et al.|[2006h . while in lO'Halloran. Satvapal fc Dudikl 
(|2006h (hereafter Paper I), we examined a sample of 18 
galaxies of varying metallicity (from 1 /50th to super so- 
lar) with high star formation rates in order to determine 
whether supernova-driven shocks do indeed play a role 
in the PAH deficit in low metallicity environments. If 
we consider the ratio of the 26 /im [Fell] line and the 
12.8 yum [Nell] line as a tracer of the strength of su- 
pernova shocks, we found a strong anti-correlation sug- 
gesting that strong supernova-driven shocks are indeed 
present within low metallicity galaxies. Furthermore, the 
PAH deficit within these objects may indeed be linked to 
the presence and strength of these shocks. However, it 
has not as yet been conclusively proved that shocks are 
the dominant process behind the PAH deficit. As already 
noted, photoionizat ion and perhap s delayed injection of 
dust into the ISM (Galliano 20Q8|) each make their own 
contributions to the PAH deficit, but as yet, it remains 
unclear as to which of these three processes is the most 
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dominant. Additionally, is the relationship between the 
[Fell] / [Nell] ratio and the PAH emission, as indicated in 
Paper I, really probing a causal effect? It could very well 
be that this relationship may be coincidental, or that 
physical or selection effects may be the primary drivers 
for the observed relationship. It would therefore be ad- 
visable to further explore the nature of the ISM within 
such star-forming environments using additional mid-IR 
probes in order to expand upon our understanding of the 
[Fell] / [Nell] ratio, and by extension, its relationship with 
key ISM indicators such as PAH strength and metallicity. 

2. OBSERVATIONS AND DATA ANALYSIS 

To accomplish this, we have expanded our sample from 
the 18 objects presented in Paper I by including archival 
IRS observations of 23 additional objects, bringing the 
sample total to 41. The full list of targets is given in Ta- 
ble 1, and listed by increasing metallicity. These galaxies 
range in metallicity from extremely low (such as I Zw 18, 
with Z I Zq— 1/50) to super-solar metallicity (> 1 Zq) 
galaxies such as NGC 7714. To differentiate between low 
and high-metallicity galaxies, we use a metallicity (12 
+ log [O/H]) cut-off value of 8.85 - just less than solar, 
as per Paper I. The galaxies range widely in morphology 
from blue compact dwarfs such as I Zw 18 to spirals such 
as NGC 7714 and IC 342. None of the galaxies in our 
sample are known to harbour AGNs. This is important, 
as PAH destruction can occur close to an AGN due to 
the hard ionization environment (e.g. Sturm et al. 2000). 
In addition [Fell] emission can be elevated in galaxies 
harbouring AGN. One possible exception is NGC 7714, 
which is optically classified as a LINER (Tho mas et al.l 
2002). However, the lack of [NeV] emission at 14 and 
24 /im and the absence of any evide nce for an obscured 
AGN by recent C handra imaging (jBrandl et all 12004 
ISmith et all 120051 ) strongly suggests that NGC 7714 is 
a pure starburst and it is therefore included in our sam- 
ple. 
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We extracted low and high resolution archival spec- 
tral data from the Short-Low (SL) (5.2 - 14.5 /im), 
Short-High (SH) (9.9 - 19.6 /im) and Long-High (LH) 
(18.7 - 37.2 /im) modules of the Spitzer Infrared Spectro- 
graph (IRS). The datasets were derived from a number 
of Spitzer Legacy, GO and GTO programs released to 
the Spitzer Data Archive, and consisted of either spec- 
tral mapping or staring observations. We obtained fluxes 
for the nuclear positions only from the mapping obser- 
vations. All the staring observations were centered on 
the galaxy's nucleus. The data were preprocessed by 
the Spitzer Science Center (SSC) data reduction pipeline 
version 15. 3 4 before being downloaded. Further process- 
ing was done within the IRS BCD-level data reduction 
package SPICE, v 2.0. 1 5 . Both the low and high reso- 
lution spectra were extracted by SPICE using the full- 
aperture extraction method. As the [SIII] emission is 
extended beyond the extent of the high resolution slits 
in particular, we applied the ALCF and SLCF extended 
source corrections to the spectra as part of the SPICE 
extraction 6 . The slit for the SH and LH modules is too 
small for background subtraction to take place and sep- 
arate SH or LH background observations do not exist for 
any of the galaxies in this sample. For the SL module, 
background subtraction was done using either a desig- 
nated background file when available or by subtracting 
the data from one nod position from the other along the 
slit. In some cases, the slit was enveloped in the host 
galaxy and background subtraction could not take place. 
For both high and low resolution spectra, the ends of 
each order where the noise increases significantly were 
manually clipped, as were hot pixels. 

For the high resolution observations, we required 
matched extractions in terms of angular extent from both 
the SH & LH slits, in order to accurately derive densities 
based upon the 18 and 33 /im [SIII] lines. The ioniza- 
tion potential of [SIII] is ~35 eV, with the [SIII] emission 
arising from gas ionized by young stars for our sample. 
With a relatively low ionization potential, [SIII] emission 
may be quite easily extended in spatial extent beyond the 
nuclear single-slit pointings for SH and LH. I ndeed, com- 
parisons of ISO and IRS data (|Dudikll2007l ) show [SIII] 
emission within star-forming AGN that extends beyond 
the size of both the SH and LH slits. Given that [SIII] 
emission is extended and we would expect to see a higher 
flux in the larger LH slit, any derived [SIII] line ratio will 
be artificially depressed. To avoid such aperture effects 
between the SH and LH slits, we scaled the LH flux by 
multiplying it by the ratio of the SH/LH slit angular 
sizes. 

As a number of the sample starbursts at high metallic- 
ity are known to be dusty systems with large obscuration, 
the extracted line fluxes were then corrected for extinc- 
tion for the sample as a whole. The fluxes and statistical 
errors, plus the derived line ratios, are presented in Ta- 
bles 1 and 2 respectively, obtained using the IDL-base d 
analysis package SMART, v. 6. 2. 6 (jHigdon et alJ 12004( 1. 
In all cases detections were defined when the line flux was 
at least 3a. For the SL, SH and LH modules the spectral 

4 Spitzer Observers Man ual, URL: 
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a URL: http://ssc.spitzer.caltech.edu/postbcd/spice.html 
6 URL: http://ssc.spitzer.caltech.edu/postbcd/doc/spice.pdf 



resolution was X/8X 60-127 and ^600 respectively, with 
FWHM in the order of 1.9 x 10" 1 /im for the PAH fea- 
ture in the low-res spectrum, and 2-4 x 10~ 2 /im for the 
fine structure lines in the high resolution spectra. For the 
absolute photometric flux uncertainty we conservatively 
adopt 15%, based on the assessed values given by the 
Spitzer Science Center (SSC) over the lifetime of the mis- 
sion. This error is calculated from multiple observations 
of various standard stars throughout the Spitzer mission 
by the SSC. The dominant component of the total error 
arises from the uncertainty at mid-IR wavelengths in the 
stellar models used in calibration and is systematic rather 
than Gaussian in nature. The 6.2 /im PAH strength 
was determined according to t he prescription adopted 
in a number of previous stud ies (Rig oupolou et alj|1999t 
iForster Schreiber et alJl2004( l. While the 7.7 and 8.6 /im 
PAH lines are stronger, the 8.6 /im line can suffer extinc- 
tion effects due to a potential silicate absorption feature 
at 9 /im, and both of these pose difficulties in determin- 
ing the continuum level. For this work, the continuum 
was determined at the center of the 6.2 /im feature, using 
a first order linear fit to the 5.5 and 11.5 /im bandpass. 
The PAH strength was then calculated as the ratio of the 
6.2 /im feature peak intensity to the underlying contin- 
uum. We assume a Gaussian fit while fitting the PAH 
and fine structure lines with SMART. 

3. RESULTS 

In Fig 1, we plot [Fell] / [Nell] vs the 6.2 /im PAH peak 
to continuum ratio to check if the relationship between 
the [Fell] / [Nell] ratio and the PAH peak to continuum 
strength first noted in Paper I, using the 6.2 /im PAH 
feature instead of the 7.7 /im, as one can perform con- 
tinuum fitting with a higher degree of confidence with 
the 6.2 /im feature, still holds for the expanded sample, 
as noted above. Employing a Sp earman rank correlation 
analysis ((Kendall fc Stuartl ll976f) to assess the statistical 
significance of this trend yields a correlation coefficient 
of r s of -0.705 and P s of 2.69 x 10 -6 , confirming the sig- 
nificant anti-correlation between PAH strength and the 
[Fell] / [Nell] ratio as seen in Paper I. The Spearman rank 
correlation technique has the advantage of being non- 
parametric, robust to outliers and does not presuppose 
a linear relation. We plot the [Fell] /[Nell] ratio versus 
the metallicity in Fig. 2, and we again see a strong anti- 
correlation (r s of -0.777 and P s of 9.36 x 10 -7 ) between 
[Fell] / [Nell] and metallicity as seen in paper I. 

Based upon what we have already seen from Paper I 
and from Figs. 1 and 2, if the [FeII]/[NeII] ratio is truly 
indicative of the strength of supernova-driven shocks 
within the extent of the high resolution slits, one would 
expect the passage of such shocks to affect conditions 
within the local ISM in quite a substantive manner. 
Shoc ks, in addition to removing dust and P AH from the 
ISM ((Reach fc Rhdl200fl iReach et al.ll2006( l. should also 
be adept in removing gas - one would therefore expect 
the propagation of intense SNe-driven shocks into the 
ISM of star forming regions to greatly impact on the 
density of the gas. In order to determine how the gas 
density within these nuclear star forming regions corre- 
sponds with the strength of the supernova-driven shocks, 
we require a reliable mid-IR tracer to probe the gas den- 
sity. The [SIII] 18.7/33.4 /im line ratio provides such a 
reliable mid-IR tracer of the gas density, as it is ideal 
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PAH 6.2 jum / 6.2 /im continuuum 

Fig. 1.— Plot of the [FcII]/[NcII] ratio versus the PAH 6.2 micron 
peak to continuum ratio for the extended sample. The observed 
trend confirms the relationships between the [Fell] /[Nell] ratio and 
the PAH strength seen in Paper I. 



[Fell]/[Nell] 26 /im / 12.8 /im 

Fig. 3.— Plot of the 18.7/33.4 fim [SIII] ratio as a function of the 
[FeII]/[NeII] ratio for the extended sample. The [SIII] ratio corre- 
lates strongly with the [Fell] /[Nell] ratio, with high values of the 
[FeII]/[NeII] ratio corresponding to low [SIII] ratios and metallicity. 
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Fig. 2. — Plot of the [Fell] / [Nell] ratio versus the metallicity 
for the extended sample. The observed trend confirms the strong 
relationship between the [Fell] / [Nell] ratio and metallicity seen in 
Paper I. 



Fig. 4.— Plot of the 18.7/33.4 fim [SIII] ratio as a function of 
the 6.2 (im PAH peak to continuum strength for the extended sam- 
ple. The [SIII] ratio correlates quite well with the PAH strength, 
indicating that the gas density scales with the increased strength 
of PAH emission. 



for probing gas (with high critical densities) in the re- 
gions surrounding starbur sts, especially as the [SIII] lines 
are starburst dominated (|Verma et alj [2004) . This ra- 
tio is sensitive to changes in the density for the 50 < 
n e < 10 4 cm" 3 , but is insensitiv e to changes in tem- 
perature (jRigopoulou et all Il996f) . In Fig. 3, we plot 
the logarithm of the ratio of the [SIII] fluxes versus the 
[Fell] / [Nell] flux ratio. There is a strikingly strong trend 
between the two line ratios, with high [Fell] /[Nell] values 
corresponding to lower [SIII] ratios. Using the Spearman 
rank correlation analysis, we get (r s ) of -0.891 between 
the [SIII] ratio and the [Fell] / [Nell] with a probability of 
chance correlation (P s ) of 7.43 x 10~ 8 , indicating a sig- 
nificant anti-correlation. Interestingly, we see a general 
decrease in the metallicity of the object with lower [SIII] 
ratios, corresponding to higher [Fell] /[Nell]. 

We again plot the logarithm of the ratio of the [SIII] 
fluxes in Fig. 4, but this time against the 6.2 micron 
PAH peak/continuum ratio. As with Fig. 3, we see a 
strong trend, but this time with the PAH strength in- 



creasing with log [SIII] - the PAH strength is increasing 
with higher gas density. Again running a Spearman rank 
correlation analysis, we get i s of 0.750 and P s of 9.64 
x 10~ 7 , confirming a significant correlation between log 
[SIII] and PAH strength. We again see that a general 
trend exists with PAH strength increasing with metallic- 
ity and increasing values of log [SIII] . 

4. DISCUSSION 

4.1. Can we explain the observed relationships? 

4.1.1. Are there other factors which may be influencing 
the observed relationships? 

As we have seen from Figs. 1-4, a series of relationships 
would seem to exist between the [Fell] / [Nell] ratio, the 
metallicity, the PAH strength and the [SIII] ratio, with 
low values of the [SIII] ratio corresponding with high 
[Fell] / [Nell] ratios and weak/no PAH emission. How- 
ever, before we press ahead with more detailed investi- 
gations of the nature of these relationships, it is prudent 
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to first check whether other factors, such as extinction 
and aperture effects might be responsible. 

4.1.2. Extinction 

As we have noted previously in Section 2, we cor- 
rected the IRS fluxes for extinction, with Ay ranging 
from ~ 2 to 30 mag. To determine the correction, we 
used the ratios of H/3/Ha and Br/3/Bra where avail- 
able from the lite r ature . Emissivity coefficients from 
IStorey fe Hummer] (|1995T ) for case B recombination with 
an electron temperature of 10 4 K and density of 100 cm -3 
were used to calculate the intrinsic line ratios for the 
H recombination lines. We note that while the den- 
sity used in these calculations is lower than that seen 
in some of our sample, we noted little change in the in- 
tr insic line ratio s for denser environments. We adopted 
a iDraind (Il989f) extinct i on law above 8.5 fim, and the 
iLutz. Veilleux fe Genzell (|1999f ) galactic extinction law 
below this point. As the observed relationships remain 
post correction, we conclude that extinction can be ruled 
out as a driver for the observed correlations. For the op- 
tical [SII] lines used lines used as comparisons to the 
MAPPINGS output, we obtained and used from the lit- 
erature extinction corrected fluxes where available. 

4.1.3. Aperture effects 

As noted earlier, the ionization potential of [SIII] is 
^35 eV, with the [SIII] emission arising from gas ionized 
by young stars. With a relatively low ionization poten- 
tial, [SIII] emission may be quite easily extended in spa- 
tial extent beyond the nuclear single-slit pointings for SH 
and LH. To compensate, we applied the SLCF and ALCF 
extended source corrections to the spectra. Additionally, 
as we are looking to derive a physical quantity, namely 
the gas density, over a similarly sized physical area given 
the extent of the [SIII] emission, it is therefore crucial 
that we use [SIII] fluxes derived from slits of equal an- 
gular size. To do this, we obtained a fraction of the LH 
full slit flux by multiplying it by the SH/LH angular size 
ratio, as noted in Section 2. The final, corrected fluxes 
are given in Table 1. 

4.2. Modeling the observed relationships 

Since we can rule out extinction and aperture effects 
as drivers for the observed correlations, we can be con- 
fident that the correlations seen in Figs. 1-4 are indeed 
true physical relationships between the strength of the 
[Fell] / [Nell] ratio and a number of ISM indicators, with 
the primary culprit for the observed relationships being 
the passage of SNe-driven shocks. However, while we 
have focussed up to now with a supernova-derived shock 
origin for the behaviour of the [Fell] / [Nell] ratio, it may 
not be the only explanation. For example, llzotov et al.l 
(2006) note sign of strong depletion of iron onto dust 
grains, and gradual destruction of those grains on a time 
scale of a few Myr, based on a survey of metal-poor galax- 
ies from the 3rd release of the SDSS. Such a process could 
undoubtedly drastically affect the nature of Fe emission, 
and by extension the behaviour of the [Fell] /[Nell] ra- 
tio, within our sample - we may instead be probing an 
abundance-driven relationship and the other observed re- 
lationships presented here would be purely coincidental 
to the [Fell] /[Nell] ratio. In order to explore if this 



indeed is the true cause or if shocks alone can explain 
the observed relationships, we used shock and standard 
Hn - PDR models in an effort to model the observed 
[Fell] / [Nell] relationship for a wide variety of environ- 
ments. We used b oth the MAPPINGS III p hotoion- 
ization/shock code ([Sutherland fe Dopitalll993h a nd the 
Cloudy photoionization code (|Ferland et al.lll998D m or- 
der to determine the relative proportions to the [Fell] 
emission from shocks and Hn/PDRs, and by extension, 
the driving process behind the behaviour of the observed 
[Fell] / [Nell] line ratio. In addition to the mid-IR out- 
put, we also wished to explore if a similar relationship 
is present at optical wavelengths. To this end, we use 
the ratio of extinction corrected [SII] fluxes at 6716 and 
6731 A taken from the literature and compared the ob- 
served ratios with the MAPPINGS output. 

4.2.1. MAPPINGS models 

In an effort to quantify the true nature of the 
[Fell] / [Nell] ratio, we initially constructed a number 
of MAPPINGS models where density, metallicity and 
luminosity were allowed to vary. The theoretical H II 
region models were generated by the MAPPINGS III 
code, which uses as input the EUV fields predicted by 
the stellar popu l ation synthesis models STARBURST99 
(jLeitherer et al.l Il999f ) . The photoionization modeling 
and shock models carried out with MAPPINGS III 
for th is anal ysis are descr i bed in iDopita fe Sutherland! 
(|1996l ) and iKewlev et all (|2001h . and are described 
briefly here. To model the Hn region spectrum, MAP- 
PINGS assumes the metallicity and the shape of the 
EUV spectrum are defined (through the STARBURST99 
input), and characterizes the local ionization state by a 
local ionization parameter: 



where S is the ionizing photon flux through a unit area, 
n is the local number density of hydrogen atoms and 
q is the maximum velocity of an ionizati on front that 
can be driven by the lo cal radiation field (jDopital [20001 : 
iKewlev fe Dopital [2002) . We used a 5 Myr continuous 
star- formation SED from STARBURST99 as input to 
the MAPPINGS III code. The parameters of the SED 
consisted of a Saltpeter IMF with a power law of 2.35 
and a star-formation rate of 1 Mq yr _1 . 

We then constructed a series of plane parallel, iso- 
baric models with P/k — 10 5 cm~ 3 K, where either the 
gas density, metallicity or the luminosity for the STAR- 
BURST99 input were allowed to vary, with the other two 
variable remaining constant. For the STARBURST99 
input models, the density, metallicity and luminosity 
were varied as follows: 

• Density: 10-1000 cm -3 ; metallicity (0.4Z Q ) and 
luminosity (10 ergs) remain constant; 

• Metallicity: 0.05 - 2Z Q ; density (100 cm" 3 ) and 
luminosity (10 ergs) remain constant; 

• Luminosity: 10 41 - 10 45 ergs; density (100 cm~ 3 ) 
and metallicity (OAZq) remain constant. 

Elemental abundances for 0.5-2.0Z Q were adopted for 
the variable metallicity models, while for the variable 
density and luminosity models, OAZq metallicity abun- 
dances were used. 
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Fig. 5.— Plots of the [SIII] (top) and [SII] (bottom) ratios vs 
[FeII]/[NeII] for the extended galaxy sample, overlaid with model 
data from MAPPINGS. The dot-dash line indicates the low density 
limit. 



Dust physics is treated explicitly through absorp- 
tion, grain charging and photoelectric heating and 
the use of a standard MRN grain size distribution 
for both carbonaceous and silic aceous grain types 
(jMathis. Rumpl fc Nordsieckl [l977h . The MAPPINGS 
c ode took the undepleted so lar abundances to be those 
of lAnders fc Grevessd ([1989) . For non-solar metallicities, 
MAPPINGS assumes that both the dust model and the 
depletion factors are unchanged. All elements except ni- 
trogen and helium are taken to be primary nucleosynthe- 
sis elements. As for the shock models, we use the pre-run 
shock models from the MAPPINGS website, which are 
based upon a grid of n = 1.0 models with v = 100-1000 
km/s in steps of 25 km/s. The magnetic parameter B = 
y/^/n is fixed to 3 in all calculations. This value corre- 
sponds to equipartition between ther mal and magnetic 
pressures (jDopita fc Sutherlandl[l996l) . 

From the MAPPINGS output from the variable metal- 
licity model, we saw that whilst the [SIII] and [SII] values 
remain relatively constant (and by extension, the den- 
sity), the [Fell] /[Nell] ratio varies quite considerably. 
With the variable density model, interestingly, we saw 
that [Fell] /[Nell] does not drop with increasing den- 
sity contrary to what we observe (Fig. 3), but rather 



slightly increases, contrary to what we observe. Finally, 
for the variable luminosity model, both the [Fell] /[Nell] 
and unsurprisingly the [SIII] and [SII] values, remained 
constant. From the MAPPINGS model data, it would 
seem that the change in metallicity is primarily respon- 
sible for the observed drop in the [Fell] /[Nell] ratio. In 
order to directly compare the photoionization and shock 
model output with the IRS and optical data, we then set 
about defining models that would replicate the sort of 
density, metallicity and luminosity parameters typically 
seen in the sample. Since the low metallicity objects tend 
to have lower densities and luminosities, we generated 
MAPPINGS models that gradually increased the metal- 
licity, density and luminosity, over the ranges outlined 
above, moving from low to high metallicity (0.05, 0.2, 
0.4 and solar metallicity), low to high density (20-1000 
cm~ 3 ) and low to high luminosity (10 42 to 10 ergs). 
As we increased the metallicity, we concurrently moved 
to high densities and higher luminosities. As with the 
earlier MAPPINGS models, we used a 5 Myr continu- 
ous star-formation SED fro m STARBURST99 as input 
to the MAPPINGS III code (jLeitherer et aLl ^g). The 
parameters of the SED consisted of a Saltpeter IMF with 
a power law of 2.35 and a star-formation rate of 1 M Q 
yr _1 . Elemental abundances for 1/20Z©, 2/5Zq and IZq 
were again adopted. 

Using the observed [SII] 6716 & 6731 A, [SIII] 18 & 
33 /mi plus the [Fell] 26 /im and [Nell] 12.8 fluxes, 
we obtained the respective ratios and performed best fit 
line calculations to the [SIII] vs [Fell] /[Nell] and [SII] 
vs [Fell] / [Nell] model data, where the line fluxes are the 
product of both the photoionization and shock model 
outputs. Using the best fit models, we plotted the model 
ratio data in conjunction with the observed IRS and op- 
tical data ratios (Fig. 5), with the model data derived 
from the second MAPPINGS run describing very well 
the behaviour of both sets of observed line ratios. 

4.2.2. Cloudy models 

The slight increase in the [Fell] emission with density 
seen from the MAPPINGS output poses the following 
question - where does this increase originate from? Is 
this increase due to an increased contribution to the 
[Fell] emission by dense Hn regions and PDRs? If 
so, what is the level of this contribution? To resolve 
this, we used the spectral synthesis code Cloudy to try 
and determine how much of a contribution H n regions 
and PDRs make to the overall [Fell] emission, how this 
contribution is affected by the gas metallicity (as sug- 
gested by MAPPINGS), and the effect of the Hn/PDR 
contribution has on the overall [Fell] /[Nell] ratio. We 
used the developmental version o f the spectral synthe - 
sis code Cl o udy, l ast de scribed bvlFerland et al.l (|1998f) . 
lAbel et all (|2005h and IShaw et alj (|2005h describe re- 
cent advances in its t reatment of PDR chemistry, while 
Ivan Hoof et al.l (|2004h describes the grain physics. As 
with the MAPPINGS models, we used a 5 Myr contin- 
uous star-formation SED from STARBURST99 as in- 
put to Cloudy. The parameters of the SED consisted of 
a Saltpeter IMF with a power law of 2.35 and a star- 
formation rate of 1 M Q yr _1 . 

Our model geometry consists of a plane-parallel slab il- 
luminated on one side by a source of UV radiation. Our 
calculations start at the hot, illuminated face of the slab, 
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where all the hydrogen is ionized. We end our calcula- 
tions at two locations, when the hydrogen ionization frac- 
tion falls below 1% (for the pure Hn region calculations, 
see below) and when the fraction of hydrogen in the form 
of H 2 ((2x n(H 2 ))/n(H) total ) exceeds 90%. Beyond this 
depth, we expect there to be little contribution to either 
[Fell] or [Nell] emission. 

We parameterize the ionizing continuum with the 'ion- 
ization parameter', which is the dimensionless ratio of 
hydrogen ionizing flux to density, 

U= Qh 
AnR 2 N H c 

where Qh is the total number of hydrogen-ionizing pho- 
tons emitted by the central object per second. For all 
calculations, we chose Logft/j = -2 5, a v alue typical in 
H II regions (jVeilleux fc OsterbrocldflQSl . The ioniza- 
tion parameter in MAPPINGS is related to the one used 
in Cloudy by the equation: 



[Fe H) 25.99nm / * [Ne n]12.81|jm 



u 



(2) 



We also include cosmic rays in our calculations. Pri- 
mary and secondary cosmic ray ioniz ation proce s ses ar e 
treated as described in Appendix C of lAbel et al.l (|2005f ). 
with an assumed cosmic ray ionization rate of 5 x 10 7 
s _1 . Cosmic rays are a significant heating source and 
also drive the ion-molecule chemistry deep in the molec- 
ular cloud. However, since our calculations stop short of 
the molecular cloud, for our models cosmic rays are only 
of secondary importance. Our calculations include the 
major ionization processes that can affect the ionization 
structure. This includes all stages of ionization for the 
lightest 30 elements. 

Dust is known to play an impo rtant role in both H n 
regions and PDRs (jDraind l2003f ). We self-consistently 
determine the grain temperature and charge as a func- 
tion of grain size and material, for the local physical con- 
ditions and radiation field. This determines the grain 
photoelectric heating of the gas, an important gas heat- 
ing process, as well as collisional energy exchange be- 
tween the gas and dust. We also treat stochastic heatin g 
of grains as outlined in iGuhathakurta fe Draind i|1989f ) , 
which can affect the dust continuum shape. We include 
grain charge transfer as a general ionization - recombina- 
tion p rocess, as described in Appendix B of lAbel et al.l 
(2005). The rates at which H2 forms on grain sur- 
faces is derived using t he temperature and ma terial- 
dependent rates given in ICazaux fe Tielensl (|2002ft . The 
assumed grain size distribution is representative of the 
star-forming regions. The ratio of total to selective 
extinction, Ry = Ay /{As- Ay), is a good indicato r 
of the size distribution of grains ijCardelri et alJll98Hh . 
ICalzetti et al.l (|2000D derived an average value for Ry 
~ 4.3. We the refore use the Ry = 4 grain size distri- 
bution given in IWeingartner fe Drainel (|2001). We also 
include size-resolved PAHs in o ur calculations, with th e 
same size distribution used by iBakes fe Tielensl (|1994f ) . 
The abundance of carbon atoms in PAHs that we use, 
nc(PAH)/ra#, is 3 x 10~ 6 . PAHs are thought to be 
destroyed by hydrogen ionizing radiation and co agulate 
in mo lecular environments (see, for instance, lOmontl 
(1986)). We model this effect by scaling the PAH abun- 
dance by the ratio of H°/H tot (n c (PAH)/n H = 3 x 10~ 6 
HH°)MH tot )]). 
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Fig. 6. — Cloudy photoionization model output of the 
[Fell] /[Nell] ratio as a function of density and Z, with contributions 
to the[FeII] emission by Hn regions only. 



The assumed gas-phase and dust abundances are 
extremely important in our calculations. We also 
assume gas-phase abundances based on an average 
of the abundanc es in th e Orion N eb ula derived 
bv iBaldwin et all (|199lD . iRubin et all l| 19911) . and 
lOsterbrock et al. (1992). The abundances by number are 
He/H = 0.095, C/H= 3 x 10" 4 ; 0/H= 4 x 10~ 4 , N/H= 
7 x 10~ 5 , and Ar/H= 3 x 10~ 6 . We also assume default 
grain abundance of Ay/iV(H tot ) = 5 x 10 -22 mag cm 2 . 
We then scale the metal and dust abundance such that 
the final value of Z = 12+log[0/ff] ranges from 7 to 10, 
in increments of 1 dex. 

There is one major exception to our assumed abun- 
dance s, and that is the Fe/H ratio. iRodriguez fe Rubinl 
( 2005) discuss the variation in Fe/H with Z, and as noted 
earlier. Hzotov et ail (|2006f ) also determined the variation 
in elemental abundances with Z from the Sloan Digital 
Sky Survey, i ncludi ng Ne and Fe. For these two elements, 
llzotov et al.l (|2006fl provide the following equations to 
derive their respective abundances (for t his work, we 
use th e notation Z instead of X, as used in llzotov et al.l 
(2006]): 



Log{Ne/0) = 0.097Z- 1.542 



and 



Log(Fe/0) = -0.606Z+ 2.994 
where Z = 12 + Log(0/ff), we find that: 
'Ne' 



Log 



Fe 



= -4.536 + 0.703Z. 



(3) 



(4) 



(5) 
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[Fe II] 25.99pm / *[Ne II]12.81jim 



0.02- 




[FeII]/[NeII] ratio as a function of density and Z with contribu- 
tions to the [Fell] emission from both H II and PDR regions. 



In all calculations, we scale the Fe abundance so that 
equation 3 is always obeyed. 

The relationship between density, temperature, and 
pressure is determined by the assumed equation of state. 
We perform two sets of calculations, one where we stop 
the calculation at the hydrogen ionization front and one 
where we include the PDR. For the pure H n region mod- 
els we assume constant density, with the density ranging 
from Logfnjj] 0.5 to Log[n#] = 4, in increments of 0.5 
dex. For the Hn + PDR calculations, we assume con- 
stant pressure, with thermal and radiation pressure being 
the dominant contributors to the total pressure. We do 
not include turbulent pressure in the equation of state, 
although we assume a small turbulent broadening of 1 
km s" 1 for its effects on line optical depths. Such an 
equation of state neglects shocks or other time-dependent 
effects. In the constant pressure scenario, we specify uh 
at the illuminated face, then solve for density based how 
the pressure contributions change with depth. Our den- 
sities range over Log[n#] = 2 to 4, in increments of 0.5 
dex. Overall, our combination of and Z represent 52 
calculations. 

In Fig. 6, we plot the [Fell] /[Nell] line ratio as a func- 
tion of Z and density for the H II region only. In Fig. 
7, we again plot the [Fell] /[Nell] line ratio as a function 
of Z and density, but this time including the PDR in 
addition to the H n region. As one can see, the relative 
H II contribution dominates over the PDR contribution 
at low densities (and metallicities), with the PDR contri- 
bution becoming more dominant as one moves to higher 
densities (~ 10 3 cm" 3 ) and metallic ities. Previous IR 
studies of starbu rst H n/PDR regions dCarral et al.lll99l 
lLord et al"lll996[ ) support this scenario. For very dense 



H II regions where the PDRs are irradiated by intense 
FUV, and thus PDR-derived [Fell] emission is dominant. 
Correspondingly at lower densities (~10 2 cm" 3 ), the Hn 
regions are larger, the PDRs lie further from the stars, 
and the resultant lowered FUV fluxes and densities do 
not excite [Fell] in the PDR, and thus the Hn region 
dominates the Hn/PDR contribution to [Fell] produc- 
tion for such environments. 

4.2.3. What is driving the Fe/Ne relationship? 

From the MAPPINGS and Cloudy simulations and the 
comparisons with the IRS and optical data (Fig. 5), 
it would seem that metallicity is the primary driver for 
the behaviour of the [Fell] /[Nell] ratio, as metallicity 
determines the origin (shocks vs Hn and/or PDRs) and 
strength of the [Fell] emission. We can reproduce the 
observed [Fell] / [Nell] dependence on Z if we assume the 
Izotov variation of Fe/Ne with Z. Cloudy reproduces the 
[Fell] / [Nell] ratio decrease with increases in density (for 
the pure H II region only) because of the following: the 
critical density for electron excitation of [Fe II] is 1.2 x 
10 4 cm" 3 , while for [Ne II] it is 5 x 10 5 cm -3 . The ratio 
of [Fe II] /[Nell] thus depends on the ratio: 

(1 + n critical (neon)/n(H}) 
(1 + n cri ticai{iron)/n(H)) 

In the low density limit, this ratio is ~42. However, for 
a density of 1 x 10 4 cm -3 , which is the high end of our 
parameter space, this ratio is only ^20 because [Fe II] 
emission is suppressed by the higher density. The reason 
the ratio increases with density in Fig. 7 (H 11+ PDR) is 
because [Fe II] is excited by atomic hydrogen in the PDR 
in addition to electrons in the H n region. The critical 
density of [Fc II] with H atoms is ~1 x 10 6 cm" 3 . There- 
fore [Fe II] emits readily in the PDR, while [Ne II] does 
not emit because all the neon is neutral. However, we 
do note that the H II and PDR contribution to the [Fell] 
emission is surpassed by the contribution from shocks 
at the lowest metallicities and densities, based on com- 
parison with the MAPPINGS output and the observed 
data - perhaps only ~10% of the [Fell] emission comes 
from H II and PDRs with that contribution rising to ~80- 
90% as one moves to higher metallicity. As we move to 
roughly solar metallicity, it would seem likely that the 
dominant mechanism for production of [Fell] (and thus 
the behaviour of the [Fell] / [Nell] ratio at high metallic- 
ities) is from H II and PDRs in combination rather than 
from SNe shocks. 

4.2.4. Can we tie in the behaviour of the gas density 
with the [Fell] /[Nell] ratio, metallicity and PAH 
strength? 

While abundances may play a role in the overall be- 
havior of the [Fell] /[Nell] ratio, the shock scenario would 
seem to fit the observed variation of the gas density and 
the PAH ratio with metallicity. Assuming that shocks 
are indeed behind the observed correlations, the follow- 
ing scenarios would seem to account for the relationship 
between the [Fell] /[Nell] and [SIII] ratios and the PAH 
strength. In very low metallicity dwarf starburst environ- 
ments, such as I Zw 18, the observed gas density is phe- 
nomenally low, in the order of ~20-100 cm" 3 . From Fig. 
3 and the models, it would seem that low density cor- 
relates strongly with the presence of strong supernovae 
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shocks, a harsh radiation environment which is lacking 
in populations of dust and PAH, as noted by the lack of 
PAH emission (Figs. 1 and 4) and generally flat slopes of 
the dust continua. It would seem likely therefore that we 
are probing the tenuous gas within a supcrbubble cavity 
that has been carved by the cumulative effects of numer- 
ous supernova-driven shocks. Previous surveys of low 
metallicity dw arf galaxies have found evidence for such 
superbubbles (|Strickland fc Steveni 2000; [Calzetti et all 
120041 : lOtt et al.ll2005l ). with ionized gas expanding at a 
rate of several tens of km/s driven by mec hanical energy 
from starbursts deposited into the ISM (jCalzetti et alJ 
l200l . 

For high metallicity environments, the picture becomes 
somewhat more complicated. Such objects are of course 
still the sites of high levels of massive star formation - the 
high metallicity starburst M82 is after all the template 
for starburst activity - but we see a much more benign 
environment for the survivability of dust and PAH. A 
number of new factors play a role in affecting star for- 
mation and its effects on the local environment. First 
of all, from the Cloudy and MAPPINGS models plus 
the observed data, we note that the dominant contribu- 
tion to the [Fell] flux no longer comes from shocks, but 
rather from Hil regions and from PDRs in particular 
- the lack of shocks propagating into the ISM provides 
a more benign, denser environment for the build-up of 
dust and PAH. Superbubbles are present as in low metal- 
licity galaxies, however supernova-driven shocks are no 
longer the dominant mechanism driving the expansion 
but rather stellar winds from massive stars. As a con- 
sequence, less mechanical energy is available to drive 
breakout and outflow and the superbubble can remain 
enclosed, especially where the gravitational po t ential of 
the galaxy is steep (|Skillman fe Bender! 119951 : iBurkertj 
l2004llCarzetti et al.ll2004l ). Given the inability of weaker 
shocks to drive out material and clear the local ISM, it 
is not surprising that we should see a buildup of metals, 
and the survivability of dust and PAH formed locally. 

4.2.5. Are abundances driving the Fe/Ne relationship? 

Whilst the SNe-driven shock hypothesis is convincing 
in terms of the enhanced [Fell] emission and correspond- 
ing low PAH emission, it is the not the only plausible 
explanation for the observed co r rela tion. As we have 
already noted, from llzotov et all (|2006l ) it has been sug- 
gested that the observed levels of [Fell] emission may 
instead be abundance driven. Using equations 1 & 2 to 
derive how the Fe/Ne ratio varies with metallicity, we 
see that the Fe/Ne ratio decreases with increasing Z in 
a fashion very similar to the [Fell] /[Nell] vs. Z corre- 
lation, suggesting that the variation in the [Fell] /[Nell] 
emission ratio with Z may be a consequence of varying 
Fe/Ne abundance w i th Z. Indeed, the slope of Fe/Ne vs. 
Z from llzotov et all (|2006l ) is very similar to the slope of 
the best fitting line from Fig 2 (m = -0.68) - our observa- 
tional dataset points to this being a possibility, and the 
MAPPINGS and Cloudy models further lend support to 
this scenario. 

However given the similarity in te rms of behav i our o f 
our observational dataset with the llzotov et al.l (|2006l ) 
results, it begs the question - what may be causing the 
change in Fe/Ne with Zl Since n eon is not depleted 
(jJuett et al.l l2006t lYao et al.l [2006h . the observed rela- 



tionship would then be due to variations in Fej H rather 
than Ne / H. So what is causing the iron depl etion? Since 
iron i s known to be depleted onto grains (jJuett et al.1 
120061; lYao et al.|[2006D . the increased variation in Fe/H 
would suggest that there is either more grain destruction, 
or less grain formation, as one moves to lower metallici- 
ties. 

If this is correct, the reason for the observed 
[Fell] / [Nell] with Z trend is that Fe is less depleted in 
low metallicity galaxies. This would have a profound 
impact on ISM studies, as the usefulness of iron emis- 
sion as a pure indicator of shocks may be compromised. 
However, this is not conclusive, and we need more obser- 
vations to deduce for certain the physical process which 
is controlling the [Fell] emission. In order to definitively 
figure out if shocks are truly responsible for the [Fell] 
emission at low metallicity, we require high resolution 
IR spectroscopy to unambiguously determine the pres- 
ence of shocks within the velocity profiles. The upcom- 
ing availability of such high resolution spectroscopic data 
from Herschel and SOFIA will be crucial to determine 
the true nature of iron emission at low metallicity and 
we look forward to furthering this line of investigation in 
a future paper. 

5. CONCLUSIONS 

Using archival Spitzer observations of 41 starburst 
galaxies that span a wide range in metallicity, we found 
a correlation between the ratio of emission line fluxes 
of [Fell] at 26 /zm and [Nell] at 12.8 //m and the elec- 
tron gas density as traced by the 18.7/33.4 fim [Sill] 
ratio, with the [Fell] /[Nell] flux ratio decreasing with 
increasing gas density. We also find a strong correla- 
tion between the gas density and the PAH peak to con- 
tinuum strength. The correlation of the [Fell] /[Nell] 
ra tio and the PAH peak to cont i nuum strength found 
in lO'Halloran. Satvapal fc D udik ( 2 0061 ) was confirmed 
for a larger sample. Using shock and photoionization 
models, we see that metallicity is the primary driver 
for the observed behaviour of the [Fell] /[Nell] ratio. It 
may very well be that the majority of [Fell] emission 
at low metallicity may be shock-derived, whilst at high 
metallicity, the [Fell] emission is dominated by contri- 
butions from H 11 and in particular, from PDR regions, 
and that at higher metallicity shocks may not play as 
significant a role in removing gas, PAH and dust from 
the ISM, unlike in low metallicity systems. However, the 
observed [Fell] / [Nell] emission may instead be follow- 
ing a metallicity-abundance relationship, with the iron 
being less depleted in low metallicity galaxies, a result 
that would have profound implications for the use of Fe 
emission lines as unambiguous tracers of shocks. Follow 
up high resolution spectroscopic observations will be re- 
quired to determine if this is indeed the case, or if the 
detected iron emission is truly tracing emission from the 
passage of shocks. 
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Table 1. IRS flux data. 



Target 


Metallicity 


F([NeII]) 


F([SIII]) 


F([FcII]) 


F([SIII]) 




12+log(0/H) 


@ 12.8 /urn 


@ 18.7 fiin 


@ 26.0 urn 


@ 33.4 ura 






(x 10- 20 W cm- 2 ) 


(x 10- 2i > W cm- 2 ) 


(x 10- 20 W cm-' 2 ) 


(x 10- 20 W cm- 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


I Zw 18 


7.2 


0.030+0.003 


0.006+0.001 


0.051+0.015 


0.011+0.006 


SBS 0335-052 


7.3 


0.014 ± 0.003 


0.013+0.003 


0.018+0.006 


0.021+0.008 


HS 2236+1344 


7.5 


0.032 ± 0.002 


0.063+0.008 


0.056+0.016 


0.112+0.011 


SDSS J0519+0007 


7.5 


0.207+0.014 


0.048+0.014 


0.155+0.056 




SBS 1415+437 


7.6 


0.068+0.013 


0.026+0.007 


0.076+0.024 


0.045+0.021 


HS 0837+4717 


7.6 


0.121+0.010 


0.020+0.005 


0.177+0.036 


0.033+0.011 


Tol 65 


7.6 


0.024+0.013 


0.011+0.002 


0.040+0.005 


0.022+0.008 


HS 1442+4250 


7.6 


0.226+0.010 


0.024+0.004 


0.127+0.041 


0.038+0.004 


Mrk 25 


7.8 


1.074+0.156 


0.847+0.104 


0.647+0.128 


0.870+0.026 


SBS 1030+583 


7.8 


0.075+0.014 


0.167+0.011 


0.028+0.007 


0.149+0.043 


Mrk 209 


7.8 


0.068+0.017 


0.070+0.024 


0.077+0.020 


0.124+0.024 


Mrk 36 


7.9 


0.281+0.071 


0.051+0.010 


0.130+0.024 


0.085+0.006 


SBS 0917+527 


7.9 


0.228+0.058 


0.115+0.024 


0.074+0.007 


0.084+0.008 


SBS 1152+579 


7.9 


0.928+0.300 


0.071+0.037 


0.305+0.018 


0.076+0.024 


Mrk 170 


7.9 


0.053 ± 0.015 


0.033+0.002 


0.048+0.006 


0.058+0.009 


UM 448 


8.0 


1.430+0.018 


0.397+0.031 


0.442+0.084 


0.563+0.099 


SBS 0946+558 


8.0 


0.053+0.015 


0.034+0.022 


0.026+0.003 


0.034+0.006 


II Zw 40 


8.1 


0.452+0.061 


1.375+0.026 


0.115+0.031 


1.274+0.304 


Mrk 930 


8.1 


0.219+0.074 


0.170+0.019 


0.112+0.026 


0.189+0.056 


Mrk 996 


8.1 


0.297+0.022 


0.271+0.044 


0.073+0.011 


0.184+0.046 


Tol 1924-416 


8.1 


0.464+0.039 


0.247+0.050 


0.058+0.001 


0.271+0.073 


II Zw 70 


8.1 


0.143+0.012 


0.178+0.041 


0.182+0.053 


0.268+0.075 


NGC 4670 


8.2 


0.482+0.173 


0.240+0.039 


0.309+0.042 


0.321+0.029 


Mrk 450 


8.2 


0.125+0.009 


0.083+0.029 


0.050+0.015 


0.096+0.059 


Mrk 5 


8.2 


0.106+0.013 


0.076+0.020 


0.075+0.019 


0.103+0.027 


Mrk 1329 


8.2 


0.305+0.015 


0.257+0.054 


0.060+0.002 


0.196+0.036 


NGC 5253 


8.2 


0.572+0.075 


1.172+0.039 


0.434+0.100 


1.749+0.249 


IC 342 


8.3 


34.239+5.285 


9.902+0.156 


5.793+0.025 


9.872+1.771 


Haro 3 


8.3 


2.016+0.084 


1.431+0.064 


0.306+0.042 


1.243+0.090 


UM 311 


8.3 


0.126+0.010 


0.0873+0.021 


0.056+0.009 


0.101+0.042 


Mrk 1499 


8.4 


0.150+0.014 


0.319+0.008 


0.069+0.017 


0.393+0.018 


UGC 4274 


8.5 


0.710+0.030 


0.402+0.036 


0.308+0.079 


0.393+0.093 


NGC 7793 


8.7 


0.318+0.077 


0.295+0.010 


0.088+0.025 


0.178+0.069 


NGC 4194 


8.8 


8.464+0.169 


2.151+0.042 


1.491+0.336 


2.151+0.466 


NGC 253 


8.9 


199.753+1.235 


42.928+0.780 


28.433+1.154 


29.012+0.416 


He 2-10 


8.9 


24.671+1.677 


9.034+0.156 


2.300+0.091 


6.072+0.990 


NGC 7714 


8.9 


8.546+0.169 


2.665+0.016 


1.012+0.121 


1.200+0.194 


M82 


9.0 


948.875+25.765 


44.638+0.047 


66.546+2.118 


18.416+0.411 


NGC 3049 


9.0 


2.800+0.353 


1.008+0.180 


0.834+0.016 


0.780+0.104 


NGC 1482 


9.1 


100.483+4.179 


13.346+0.104 


1.759+0.073 


4.905+0.119 


NGC 2903 


9.3 


84.257+1.871 


4.456+0.026 


7.563+0.057 


1.979+0.267 



Note. — Columns: (1) Common source names; (2) Metallicity of the galaxy; (3) Extinction corrected flux of the 12.8 fim. [Nell] fine structure 
line in Watts per centimeter squared; (4) Extinction corrected flux of the 18.7 fim [Sill] fine structure line in Watts per centimeter squared; (5) 
Extinction corrected flux of the 26.0 fim [Fell] fine structure line in Watts per centimeter squared; (6) Extinction corrected flux of the 33.4 /xm 
[SIII] fine structure line in Watts per centimeter squared. 



Tracing the [Fell] /[Nell] ratio within star forming dwarf galaxies: a Spitzer IRS archival survey 

Table 2. IRS and optical flux ratios. 
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Target PAH 6.2 micron P/C [Fell] / [Nell] log QSIII] 18.7/33.4 /an) log ([SII] 6716/6731 A) 



(1) 


(2) 


(3) 


(4) 


(5) 


T 7w 1 & 

1 Zj W lo 


1 2 


1 7 A 


n 9Q 


1 32 a 


DDO UOOO-UdZ 


1 9^ 
l.ZO 


1 9Q 
i.zy 


-0 22 


1.2 b 


HS 22^64-1 '344 

1 1 1_' iiiJU~ lO'l'l 


1.88 


1.75 


-0.25 


1.15° 


SDSS J0519+0007 


1.89 


0.75 


1.17 c 




SBS 1415+437 


1.57 


1.11 


-0.24 


1.15 b 


HS 0837+4717 


1.66 


1.47 


-0.22 


l.l c 


Tol 65 


1.41 


1.68 


-0.29 


1.21 c 


HS 1442+4250 


1.81 


0.56 


-0.20 


- 


Mrk 25 


2.56 


0.60 


-0.11 


- 


SBS 1030+583 


1.67 


0.37 


0.05 


1.25 b 


Mrk 209 


2.27 


1.12 


-0.25 


1.28 a 


Mrk 36 


1.82 


0.46 


-0.23 


1.07 b 


SBS 0917+527 


2.1 


0.33 


0.04 


1.14 b 


SBS 1152+579 


2.41 


0.33 


-0.14 


1.17 b 


Mrk 170 


2.2 


0.90 


-0.25 




UM 448 


1.63 


0.31 


-0.15 


1.07 b 


SBS 0946+558 


2.5 


0.50 


-0.01 


l.ll b 


II Zw 40 


2.6 


0.26 


0.03 


1.15 b 


Mrk 930 


2.5 


0.51 


-0.04 


1.15b 


Mrk 9Q6 


3.3 


0.24 


0.17 




rp.l i HO /I /lift 

loi iyz4-4io 


9 Q 
Z.O 


U. lo 


-U.U4 


1 1 9b 
1. 10 


II Zw 70 


2.1 


1.27 


-0.18 


1.2 b 


INCtC 4670 


3.75 


0.64 


-0.13 




Mrk 4oU 


2.64 


0.40 


-0.07 


1.17 


Mrk 5 


2.91 


0.71 


-0.13 


1.11 


ivirK lozy 


O.OO 


U.4U 


U.U t 


i.uo 




3.21 


0.71 


-0.13 


0.91 b 


1L' o4Z 


Q A 9 


n i7 
U. 1 / 


n nn 
U.UU 


U.00 


xidro o 


4. 1(5 


n it; 
u. io 


ft ftfi 


1 9^ a 

1 . zo 


UM 311 


3.27 


0.44 


-0.06 


0.92 b 


Mrk 1499 


4.27 


0.46 


-0.09 




UGC 4274 


3.92 


0.43 


0.01 


0.9 e 


NGC 7793 


4.78 


0.28 


0.22 




NGC 4194 


3.66 


0.18 


0.00 




NGC 253 


3.05 


0.14 


0.17 




He 2-10 


4.1 


0.08 


0.17 


1.04 f 


NGC 7714 


4.22 


0.12 


0.35 


0.88 b 


M82 


5.07 


0.07 


0.38 


0.8 e 


NGC 3049 


5.23 


0.30 


0.11 


0.87 b 


NGC 1482 


3.7 


0.02 


0.43 




NGC 2903 


5.26 


0.09 


0.35 


1.11 s 



Note. — Columns:(l) Common source names; (2) PAH 6.2 micron peak to continuum ratio; (3) [FcII]/[NcII] line ratio; (4) Logar i thm o f the 
[Sill] 18.7 to 3 3.4 m icro n line ratio; ( 5 ) Log arith m of the [SII] 67 1 6 to 6 731 A line ratio. So urce s for [SII] ratios (a): IKong fe Chend (1199911 : (b) 
IHovos fe Dial d2006l) : (c) IIzotov et al.1 ll2004f) : fdl llzotov fe Thuanl 11200211 : fel lHo et all dl997T) : fel lHo et all J1997H : (el lvacca fe Contil lil992fl . 



